Photonic Band Gaps in One-Dimensionally Ordered Cold Atomic Vapors 
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We experimentally investigate the Bragg reflection of light at one-dimensionally ordered atomic 
structures by using cold atoms trapped in a laser standing wave. By a fine tuning of the periodicity, 
we reach the regime of multiple reflection due to the refractive index contrast between layers, yielding 
an unprecedented high reflectance efficiency of 80%. This result is explained by the occurrence of a 
photonic band gap in such systems, in accordance with previous predictions. 
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Cold atomic vapors can be used as suitable optical me- 
dia for a number of applications or fundamental studies, 
in particular, in the fields of nonlinear and quantum op- 
tics, but also as complex optical media to study exotic 
wave transport phenomena, for instance coherent multi- 
ple scattering [l[ , weak localization 0] or random lasing 
3|. In these examples, an essential property is the disor- 
der inherent to atomic vapors that are simply confined 
in a magneto-optical trap (MOT). 

On the contrary, ordered clouds of cold atoms should 
exhibit different light-transport properties. The appear- 
ance of photonic band gaps (PBGs) has indeed been 
predicted in one-dimensional lattices 0, [I| and recent! 
in three-dimensional, diamond (non-Bravais) lattices 
The realization of PBGs in cold-atom samples would 
open the way to study new regimes of light transport 
in atomic vapors, where correlations and long-range or- 
der play a dominant role. The crossover regime, between 
order and disorder, or correlated disorder, is also a rich 
subject (see, e.g., 7|), in particular, in relation to Ander- 
son localization [8|. 

However, creating efficient photonic structures is hard 
with cold atoms because these dilute systems have a 
low refractive index [Fig. Ufa)] and a limited length. 
So far, efficient Bragg reflection of light has only been 
obtained with hot vapors, using an electromagnetically 
induced grating Previous investigations with cold 
atoms trapped in optical lattices have reported Bragg 
scattering with low efficiency. A first series of experi- 
ments used three-dimensional quasiresonant lattices fiol — 
Il2j ]. with a lattice geometry that does not create PBGs 
6]. Efficiencies below 1% were reported. A second series 
of experiments investigated the ID case [l3-15|, where a 
band gap is expected, but the maximum reflectance was 
only 5% in total power (30% corrected from the partial 
overlap between the probe beam and the atomic sample) 
[loT ]. The limitations came mainly from the probing an- 
gle, which limited the interaction length with the lattice, 
and losses, i.e., out-of-axis scattering, due to the imagi- 
nary part of the atomic polarizability near resonance. 

In this Letter, we report the observation of an efficient 
Bragg reflection at a one-dimensional lattice, reaching 



the regime of multiple reflections due to the refractive 
index contrast between layers and an 80% efficiency in 
total power. This high efficiency is obtained by (1) us- 
ing a small-diameter probe beam and a small probing 
angle to optimize the overlap with the atomic grating, 
and (2) adjusting the periodicity of the lattice so that 
the Bragg condition is fulfilled off the atomic resonance, 
thus strongly reducing the losses. Our experimental ob- 
servations are explained by the appearance of a PBG, 
which we show to be robust against the system imperfec- 
tions (finite length, varying density). 

The experiment starts with a vapor-loaded MOT of 
87 Rb containing about 6 x 10 8 atoms. A dipole trap is 
generated by a homemade titanium-sapphire laser, whose 
available power is 1.3 W and whose wavelength Adi p is 
tunable. The beam is focussed on a waist (1/e 2 radius) 
Wdip = 220 /.im at the MOT position (Rayleigh length 
zr ~ 0.2 m). A ID optical lattice is made by retrore- 
flecting the beam, thus generating a structure whose pe- 
riodicity is Adip/2. 

After stages of compression and molasses, the MOT is 
switched off and a waiting time of a few ms allows the 
untrapped atoms to fall down. Then, we can characterize 
the trapped sample with absorption imaging or acquire 
transmission and reflection spectra. Typical numbers for 
the trapped atoms are N — 5 x 10 7 atoms distributed 
over a length L ~ 3 mm (~ 7700 layers). The temper- 
ature, T ~ 100 /iK, is related to the potential depth Uq 
by a constant factor r\ = Uo/k^T ~ 3.5. The transverse 
extension of the cloud is then a± = w dip / (2 ^Jrj) ~ 60 /im 
and the thickness of each layer along the lattice axis z is 
it 2 = Adip/ (2^-^/27/) — 47 nm (rms radii in the harmonic 
approximation) . 

To acquire spectra, we shine a weak and small (waist 
w = 35 /xm) , linearly polarized probe beam onto the 
lattice under an angle of incidence ~ 2°. The an- 
gle is small enough so that the probe interacts with the 
lattice over its entire length and that the reflection is 
specular 14). The transmitted and reflected beams are 
then recorded with avalanche photodiodes. The probe 
frequency uj is swept in the vicinity of the atomic reso- 
nance uio (F = 2 — > F' = 3 closed transition of the D2 
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FIG. 1. (Color online), (a) Refractive index n as a function of the normalized detuning 8/T for an averaged density p = 
7 x 10 11 cm -3 . The dashed lines represent the Bragg condition Eq. (fTJ) with AAdip = 0.24 nm. (b) Experimental reflection 
spectra for different AAdip. (c) Measured maximal reflectance as a function of AAdip for two different atom numbers, N = 2 x fO 7 
(squares) and N = 6 x 10 7 (dots). The spectra of (b) are taken from the latter, (d) Simulated reflection spectra. 



line, A = 780.24 nm, linewidth T/2n = 6.1 MHz) by 
using an acousto-optical modulator in double-pass con- 
figuration. The other hyperfine levels are far enough to 
be negligible. The presented data are the result of an av- 
erage of typically 100 cycles (the duration of each cycle 
is ~ 1 s). 

Reflection occurs in the vicinity of a Bragg condition: 
the difference between the incident probe wavevector and 
the reflected wavevector must equal the lattice vector, 
i.e. 2n(5)kocos9 = K\ &t , where ko = 2tt/\q is the probe 
wavevector in vacuum, n is the real part of the average 
refractive index of the medium, 5 = u> — ujq is the probe 
detuning from the atomic resonance, and K\ &t = 47r/Adi p 
is the lattice vector (l6| . 

Experimentally, we keep the angle constant and adjust 
the Bragg condition by tuning the wavelength of the lat- 
tice beam. It is thus meaningful to rewrite the Bragg 
condition in the following form, 

^)-l = "^, (1) 

where AAdip = Adi P — Adi P o is the shift from the "geo- 
metric" (with n = 1) Bragg condition Adi P o = Ao/cos#. 
With 6 = 2°, A dip0 ~ 780.7 nm. Then, for a given lattice 
wavelength, the Bragg condition is fulfilled for probe de- 
tunings S given by Eq. (p}, see Fig. Ufa). There are in 
general two such frequencies, but one is almost on reso- 
nance, where losses prevent any efficient reflection. The 
other Bragg frequency [5b in Fig. [Ha)] may be farther 



from resonance and can be tuned in order to search for 
an optimum. 

Such an experiment is reported in Fig. Gib), which 
shows a set of spectra for different AAdip- As expected 
from Eq. Jl}, the spectra display a strong asymmetry, 
which evolves as the lattice wavelength is changed, the 
maximal efficiency going from one side of the resonance 
to the other side while AAdip changes its sign. One 
can also clearly observe an optimum value of AAdip 
for reaching high efficiencies, namely i? max — 80% for 
AAdip — 0.24 nm with our best atom number [Fig.QJc)]. 
Note that for each AAdip, we adjust the power accord- 
ingly to keep constant the potential depth and subse- 
quently the atom number and temperature. 

The existence of an optimum can easily be understood 
by considering the limiting cases. When AAdip is very 
small, the Bragg condition is fulfilled where the refractive 
index is almost one, in virtue of Eq. ([1]), very far from 
resonance, leading to a very small index contrast in the 
periodic structure and thus to an inefficient reflection. 
In the opposite limit, when AAdip is large, the Bragg 
condition is fulfilled near S = ±r/2, where losses, due 
to the imaginary part of the atomic polarizability, are 
large and play a detrimental role. Ultimately, if AAdip 
is too large for the given averaged atomic density p, the 
Bragg condition cannot be fulfilled at all. The optimum 
AAdip depends on the average density p: a larger density 
allows us to increase AAdip, and thus the index contrast, 
without shifting the Bragg frequency towards resonance, 
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FIG. 2. (Color online). Spectra for different atom numbers N 
in the lattice (constant length, varying density) with AAdip = 
0.15 nm. Inset: Maximum reflectance as a function of N. 



i.e. without increasing losses. This is illustrated in Fig. 
[ljc), which shows that the optimum AAdip is larger with 
more atoms. 

We have also studied the reflection spectrum as a func- 
tion of the atom number for a given AAdip- In order to 
vary the atom number while keeping the lattice length 
constant, we have changed the waiting time in the lat- 
tice, before the spectra acquisition, from 5 ms to 400 ms, 
thus varying atom losses. The result is reported in Fig. 
[2 As expected from Eq. (TTJ), we observe a shift farther 
from resonance together with a broadening for increasing 
densities. The evolution of the maximum reflectance as 
a function of the atom number (inset of Fig. [2]) reveals 
that we almost reach saturation. 

With such a high efficiency, it is natural to ask whether 
a photonic band gap occurs in our system and is respon- 
sible for the high reflectance. It is well-known that peri- 
odic ID systems give rise to PBGs for any nonzero index 
modulation 17|, so that no distinction is usually made 
between a ID PBG and a Bragg reflector. However this 
is true only for infinite, lossless, and perfectly periodic 
media. As already noted, the atomic polarizability is 
complex so that our system has losses, and it is of course 
finite. Moreover, the lattice is not perfectly homogeneous 
but has a smooth density variation, so that the period- 
icity is not perfect. The situation is thus more intricate 
and deserves a careful analysis. 

We can get a first insight with orders-of-magnitude es- 
timations. A Bragg reflector made of repeated pairs of 
dielectric layers gives rise, in the stop band, to an evanes- 
cent wave whose penetration length (or Bragg length) 
is given by L cw = Ao/(4An), where Art is the refrac- 
tive index difference between the two materials fljj . 
By approximating the Gaussian atomic distribution in 
each well by a single layer of constant density with the 
same rms width and using Eq. (TT]), the index contrast 
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FIG. 3. (Color online). Modeling of our system, (a), (b) Real 
(a) and imaginary (b) parts of the effective wavevector k c g in 
the medium as a function of the probe detuning S. The axes 
are reversed so that the graphs look like standard dispersion 
relations oj(k). (c) Normalized local density of states, (d) Re- 
flection spectrum. In all plots, the light gray line corresponds 
to a homogeneous, nonperiodic medium of the same average 
density, the dashed line to an infinite lattice and the solid red 
(medium gray) line simulates the actual density distribution 
of our system. 



is An = n^rj/Q x AAdi P /Adi P - In our experiment, the 
maximum reflection is obtained for AAdip ~ 0.24 nm 
and we have r\ ~ 3.5, which gives a penetration length 
^cw ~ 0.26 mm, sensibly smaller than our optical lat- 
tice [L ~ 3 mm), so that we are indeed in the multiple- 
reflection regime. The effect of losses can also be evalu- 
ated by comparing the corresponding attenuation length 
Li oss to L cw . Given our estimated averaged density 
p r~j 7 x 10 11 cm~ 3 and the detuning 6b determined by 
the Bragg condition (JTJ, we estimate Li oss = 1/ P<t sc (6b), 
where ct sc (5b) is the scattering cross section, which gives 
^loss ~ 3.8 mm, also much larger than L cw . Thus, this 
simple calculation confirms that our system fulfills the 
conditions L cw -C L,L\ oss , which are necessary for the 
appearance of a band gap. A number of effects are how- 
ever not taken into account, in particular the actual den- 
sity distribution along the lattice, so that a more precise 
modeling is still valuable. 

We use the transfer matrix formalism to simulate the 
wave propagation in our system, and we refer to (US EH 
for detailed descriptions in similar contexts. From the 
transfer matrix of one single period, which takes into ac- 
count the Gaussian atomic distribution in each well , 
we draw the dispersion relations of the medium (effec- 
tive wavevector k c g vs oj), valid in the limit of an in- 
finitely long lattice (Bloch theorem). These are shown 
in Fig. |3Ja,b), for the typical parameters of the experi- 
ment p = 7 x 10 11 cm -3 , i] = 3.5 and AAdip = 0.24 nm. 
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The imaginary part of fc e ff is composed of one Lorentzian 
of width r centered on resonance, which corresponds to 
losses, and a supplementary part, leading to an evanes- 
cent wave, corresponding to a band gap [Fig.[3Jb)]. In the 
same frequency range, the real part of k c g has a reduced 
variation with ui [Fig. [3Ja)], corresponding to a reduced 
density of states (DOS) V = d(Re{k cS )) / duj . This last 
formula is however not always valid, especially in anoma- 
lous dispersive media, and we use the method of [l9| to 
compute the normalized local DOS in the middle of the 
lattice. By using the complex reflection coefficients r\ , r 2 
of the two surrounding, finite or infinite [3] semilattices, 
we obtain 



V = Rc 



2 + ri + r 2 
1 - rir 2 



- 1 



(2) 



The result, shown in Fig. [31(c) , exhibits clearly the band 
gap, with a strong reduction of the normalized local DOS, 
which reaches at minimum 6% for an infinite lattice. 
The corresponding reflection spectrum reaches 85% [Fig. 
|3jd)]. In both cases, the limitation (for achieving high 
reflectance or low DOS) is the remaining losses. 

To take into account the limited length of the sam- 
ple and its actual, smooth density distribution along the 
lattice axis, that we characterized by absorption imag- 
ing, we compute the reflection spectra through the whole 
structure by multiplying elementary transfer matrices 
computed with the corresponding local density. The local 
DOS then exhibits a smaller reduction, reaching at min- 
imum 23%, but the maximum reflectance is almost as 
high as in the infinite limit, reaching 82% [Fig. [3jc,d)]. 

Finally, to simulate our experimental spectra precisely, 
we must also take into account the trapping-induced in- 
homogeneous light shift and the finite transverse size. 
Indeed, we only considered so far infinitely extended lay- 
ers, which is a crude approximation, the atom cloud hav- 
ing an rms width of 60 /im. The probe beam, having 
also a finite size, probes a distribution of local density, 
which induces an inhomogeneous broadening of the spec- 
tra. There is in fact a tradeoff for the probe beam size: 
a small beam allows us to probe a well-defined and max- 
imum density, but since the reflection is very sensitive 
to the angle of incidence, the probe divergence induces 
also an inhomogeneous broadening. Experimentally, we 
have tried several probe sizes and obtained the maximum 
reflectance with wq — 35 /im. To simulate these effects 
we averaged many spectra over Gaussian distributions of 
incident angles and densities corresponding, respectively, 
to the probe divergence and atomic density distribution 
along the transverse directions. We obtain the spectra of 
Fig. QJd), in good agreement with the experiment [20j | . 

To summarize, we have studied experimentally the 
Bragg reflection of light at a one-dimensional atomic 
quasiperiodic structure in the regime of multiple reflec- 
tion, demonstrating a high efficiency, thanks to a fine 
tuning of the lattice periodicity. Then, motivated by the 



modeling of our system, we have investigated the effects 
of the finite length and of the smooth density variation 
along the lattice on the appearance of a band gap and 
its quality. These imperfections, inherent to cold-atom 
systems, do not destroy the band gap which, even if the 
DOS reduction is not dramatic, still substantially modi- 
fies the transport properties of light, as demonstrated by 
the observed high reflectance. 

Improving further the quality of the PBG requires re- 
ducing the losses by creating the band gap farther from 
resonance, which could be done by increasing the den- 
sity or the lattice length, for example with a larger ini- 
tial MOT. Another fascinating possibility is to tailor the 
atomic polarizability to remove its absorptive part, while 
enhancing the refractive index, following the ideas of 
21 , 22| . Finally, the extension to more evolved lattice ge- 
ometries, for example with bichromatic lattices [23| and 
to three dimensions @ , has been proposed recently. This 
last case is of course the most appealing since a 3D band 
gap would profoundly modify the atom-light interaction 
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